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Abstract 

O 

oo t 

ON \ The neutral Higgs boson is expected to have a mass in the region 90-150 

D \ GeV/c 2 in various schemes within the Minimal Supersymmetric extension to 

' the Standard Model. A first generation Muon Collider is uniquely suited 

to investigate the mass, width and decay modes of the Higgs boson, since 

•i-H . 

the coupling of the Higgs to muons is expected to be strong enough for it 



to be produced in the s channel mode in the muon collider. Due to the 
narrow width of the Higgs, it is necessary to measure and control the energy 
of the individual muon bunches to a precision of a few parts in a million. We 
investigate the feasibility of determining the energy scale of a muon collider 
ring with circulating muon beams of 50 GeV energy by measuring the turn 
by turn variation of the energy deposited by electrons produced by the decay 
of the muons. This variation is caused by the existence of an average initial 
polarization of the muon beam and a non-zero value of g — 2 for the muon. We 
demonstrate that it is feasible to determine the energy scale of the machine 
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with this method to a few parts per million using data collected during 1000 
turns. 



I. THE METHOD 

The spin vector S of a muon in the muon collider will precess according to the following 
equation, first derived by Bargmann, Michel and Telegdi IJ 

f = (ixS (1.1) 
dt 

ft = — I (1 + a^)B L + (1 + a)B\\ - (07 + t^—)P x - I (!-2) 

7^/! \ 1 + 7 c y 

where B± and B\\ are the transverse and parallel components of the magnetic field with 
respect to the muon's velocity (3c, e is the electric charge, m M the mass of the muon, a = 3z ^ L 
is the magnetic moment anomaly of the muon and 7 and g are the Lorentz factor and the 
gyromagnetic ratio of the muon. The value of a = for the muon is 1.165924E-3 ||. In 
what follows, we will consider the ideal planar collider ring case where B\\ and E are zero. 
For such a collider ring, Q is given by 

Q = ft cyc (l + aj) (1.3) 

where Q cyc is the angular velocity of the circulating beam. From this, it follows that when the 
beam completes one turn, the spin will rotate by a further 07 x 2tt radians. We will compute 
the precision with which 7 can be determined by measuring the energy of the electrons 
produced by muon decay in this ideal case. We will examine the effects of departures from 
the ideal case in the last section. 

It can be shown that the angular distribution of the decay electrons in the muon center 
of mass is given by the relation || 

d 2 N 



dxdcosO 



N(x 2 (3-2x) -Px\l -2x)cos6) (1.4) 
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where N denotes the number of muon decays, x = 2E/m ll is the electron energy E in the 
muon rest frame expressed as a fraction of the maximum possible energy (« 0.5m M ), cosd 
is the angle of the electron in the muon rest frame with respect to the z axis which is the 
direction of motion of the muon in the laboratory and P is the product of the muon charge 
and the z component of the muon polarization. The muon polarization is defined as the 
average of the individual muon unit spin vectors over the ensemble of muons considered. 
We note that the distribution is linear in P. 

A routine was written to generate muon decays according to equation |1.4| . Figure [I] shows 



the shape of the function in equation |1.4j and the generated events in x, cosd space for various 
values of P. There is excellent agreement between the theoretical shape of the function and 
the Monte Carlo generated events. The average energy < E > and longitudinal momentum 



< Pi > of the electron in the muon rest frame can be obtained using equation 1.4 as follows. 



m u f f d 2 N , , „ 7m„ , 

< E >= — - / / x— -dxdcosO = £ 1.5 

2 J J dxdcosO 10 2 ( J 

m u f f „ d 2 N , , „ P m,, , 

< P L >= — ^ / / xcosO— -dxdcosO = £ 1.6 

2 J J dxdcosO 10 2 v ; 

These two quantities form the components of a 4-vector, whose transverse components are 

zero, which may be transformed to the laboratory frame to yield the average electron energy 

< ElaJb >■ 



<Elab >=^E,(l + ^P) (1.7) 



where E^ is the energy of the muon beam . Since the polarization P precesses from turn to 
turn by the amount u = 7(5 — 2)/2 x 2tt radians, and the number of muons decrease turn by 
turn due to decay and losses, the total energy E(t) due to decay electrons observed during 
turn t in an electromagnetic calorimeter will have the following expression 

E (t) = Nc(- at \^E,(l + ^(Pcosu J t + <p))) (1.8) 

where N is the number of muon decays sampled in turn 0, is an arbitrary phase containing 
information on the initial direction of polarization and a is the turn by turn decay constant 
of the muon intensity which in the absence of losses other than decay is given by 
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FIG. 1. The top lego plots shows the generated events and the theoretical decay function in 
the x,cos6 plane for P = -1.0. The lego plots at the bottom of the figure show the corresponding 
plots for P = 1.0. 
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FIG. 2. (a)Total energy observed as a function of turn number for P = -1.0 with individual 
electron energies in the range 0-10 GeV for 100,000 muon decays, (b) Electron energies in the range 
10-25 GeV (c) 25-50 GeV (d) All electrons included. Superimposed is a functional form defined by 
equation [L^ 
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where t circ is the time taken to circulate around the storage ring and tuf e is the muon life 
time. 

For a 100% polarized beam, the amplitude of the oscillations is only 1/7 that of the 



non-oscillating background. It can be seen from equation |1.4| that the sensitivity to P is 
enhanced by selecting larger values of cosO. This implies selecting electrons with higher 
laboratory energy. Figures 0(a-c) show the deposited electron energy as a function of turn 
number for polarization P = 1.0 for individual electron energy ranges of 0-10 GeV, 10-25 
GeV and 25-50 GeV respectively as a function of turn number. Figure 2(b) shows very little 
oscillatory signal, since the electrons in that energy range have small values of cosO. Figure 
2(d) shows the deposited electron energy with no electron energy cuts. Superimposed is 



the predicted behavior according equation |1.8| . This serves as a consistency check for our 



routines. The signal to background ratio increases as we demand electrons with higher 
value of cosO. In what follows, we use electrons with energy greater than 25 GeV during 
the investigative phase of this analysis and will later optimize this cut. In practice, we can 
select electrons with energies above a value by momentum analyzing them with a dipole 
field before they enter the calorimeter. 

The method to determine the energy scale of the collider would then entail fitting a 
functional form of the type 

f(t) — Ae~ Bt (Ccos{D + Et) + F) (1.10) 

to the energy observed in the calorimeter. The variables A, B, C, D, E, F are parameters to 
be fitted. The information on the energy scale is contained in the parameter E. 

A. Parameters of a 50 GeV idealized muon storage ring 



In order to arrive at reasonable numbers for a and u, we consider a storage ring of 50 
GeV muons with a uniform bending field of 4 Tesla. This would produce a circular ring 



Parameter 


Value 


Parameter 


Value 


Muon Energy 
spin precession in one turn 
radius of ring 
dilated muon life time 


DU (jG V 

3.4667 radians 
41.66666 meters 
0.10397E-02 sec 


7 

Magnetic field 
beam circulation time 
turn by turn decay constant 


A 99 

4.0 Tesla 
0.87327E-06 sec 
0.8399E-03 



TABLE I. Parameters of an idealized muon storage ring 

with the parameters given in table |. It should be noted that for an idealized storage ring 
with constant B field considered here, a does not depend on 7, since 

t (111) 

a = 1.12 

0.3Bct life v ; 

where m M is the muon rest mass, B is the bending field of the storage ring and c is the 
velocity of light. A 100 GeV collider ring will have the same a as a 50 GeV collider ring or a 
25 GeV collider ring in this idealized case. As 7 changes slightly, t circ changes in proportion, 
a being the constant used to convert measurements of t C i rc to 7. Measuring the decay rate 
of muons also affords a second method to determine 7. The beam circulation time t C i rc can 
be measured to precisions of the order of a part in 10 6 and the fractional error in muon 
lifetime is 1.82E-5 0. The fractional error in 7 obtainable by observing the rate of decay 
of the muons will then be dominated by the precision that one can measure a, namely 
57/7 = 5a/ a. 



B. Generation of events and fitting for 7 



Since equation |0] is linear in P, the decay distribution of an ensemble of muons depends 
only on P, the ensemble average of the z component of the individual muon spin vectors. 
However, because of the momentum spread of the muons, each individual particle will have 



a 7 slightly different from the average and hence the precession of the spin vector around the 
ring will be different, leading to a slightly different value of P for the next turn. We model 
the beam by generating an ensemble of 100,000 muons each having its own spin vector and 
momentum. In an actual collider, it will be possible to sample significantly more decays than 
this. During each turn, we decay all the beam particles once and record the number and 
total energy deposited by electrons with individual energies above 25 GeV. Approximately 
27% of the decay electrons pass this cut, on average. We decrease the number of decays by 
the appropriate number expected by muon decay alone for the next turn. At this stage we 
do not introduce fluctuations in the number of decays from turn to turn, since the 100,000 
muons are meant to be representative of a much larger number in the actual ring. We precess 
the 100,000 spin vectors by their individual precession rates and make them decay again. 
We repeat this for 1000 turns. We re-use the muons after each turn since the 100,000 muons 
represent our model of the muon ensemble in the collider. 

1. Generation of muon spin vectors 

We generate 4 different samples of events with different ensembles of spin vectors. The z 
component of the unit spin vector of a muon S z is allowed to vary from -1 to 1. This range is 
divided into 51 bins and the z components are generated using a binomial distribution whose 
average value is specified. We are justified in treating this problem in this classical fashion, 
since each "muon" represents an ensemble of actual muons with quantized spin components. 
A more realistic generation of the spin vectors with correlations between momentum spread 
and P would require a detailed modelling of the pion decay and muon transport systems and 
is not warranted here since the effect due to the distribution in S z is expected to be small. 
Figures (a-d) show the distributions of S z for the 4 samples. The average value of the 
distributions are 0.9, 0.74, 0.5 and 0.26 respectively. We study negatively charged muons 
resulting an initial value of P of -0.9,-0.74,-0.5 and -0.26 respectively for these samples. In 
the absence of momentum spread, the decay distributions would only depend on P and not 
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on the details of the distribution of S z . The angles of the spin vectors are precessed by the 
individual 7 dependent precession rate from turn to turn. In what follows, we assume a 
beam energy spread of 0.03% for the muons for all samples unless otherwise specified. 



2. Fitting procedure and generation of errors 



The energy deposited every turn is fitted to the functional form given by equation |1.1C 
using the CERN program MINUIT |4[]. In order to study the variation of the fractional 
error £7/7 with the number of electrons sampled, we fluctuate the energy observed in the 
calorimeter E m by 

4(1-03153) (1.13) 



< E m > 2 N 

where N is the number of electrons sampled. See Appendix for a derivation of this formula. 
We analyze the case for 41261, 10315, 2579 and 1146 electrons sampled which corresponds 
to a fractional error in the measured total energy of PERR = Q f .5E-2,1.0E-2,2.0E-2 
and 3.0E-2 respectively. 



II. RESULTS 

We simulate the muon collider spin precession for a grid of values of P =-0.9,-0.74,-0.5 
and -0.26 and fractional measurement error for the first turn (PERR) of 0.5E-2, 1.0E-2, 
2.0E-2 and 3.0E-2. Figure |(a) shows the result of the MINUIT fit plotted for 50 turns for 
P=-0.26 and PERR=0.5E-2. Figure f|(b) shows the same plot but with the function being 
plotted only at integer values of the turn number t. A beat is evident in both the theoretical 
curve and the simulated measurements as a result of sampling the oscillation function at fixed 
intervals, not connected with the oscillation frequency. The origin of the beat is stroboscopic. 
Figure |](c) shows the pulls, defined as {data — fit) / error at each measurement as a function 
of turn number for 1000 turns. There are no major turn dependent variations in this quantity 
indicating that the fit converged satisfactorily. Figure ^(d) shows the histogram of the pulls, 
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FIG. 3. (a)-(d) show the distribution of the z component of the spin vectors for the four samples 
considered. 
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which approximates a unit Gaussian as desired. Table shows the results of the fit for the 
grid of values of P and PERR. The results presented in table [Tl] are shown graphically in 
Figure^. As an example, for an average polarization P = -0.26, the fractional error in £7/7 
varies from 5.1E-6 to 1.9E-5 as the fractional error in the electron energy sampled varies from 
0.5E-2 to 3.0E-2, corresponding to the number of electrons sampled during the first turn 
varying from 41261 to 1146. The average number of decays in the muon collider is expected 
to be 3.2E6 decays per meter for a beam intensity of 10 12 muons. The error in determining 
7 is thus going to be dominated by the fluctuations in the number of electrons sampled turn 
by turn, rather than sampling fluctuations in the calorimeter. We have simulated conditions 
involving « 40,000 decays. It should be possible to go to higher statistical precision than 
computed here by sampling larger number of electrons. 

The results for £7/7 obtained from the measurement of the turn by turn rate of decay 
of the electron energy are not competetive with the precession method primarily because 
of the small value of a (0.8399E-3). This leads to larger fractional errors for 7 from this 
method (which also assumes that the loss of intensity is entirely due to the decay process) 
by almost three orders of magnitude than from the precession method. 

A. Variation of ^7/7 as a function of muon energy 

The spin precession per turn equals 27r for a 7 value of 857.689, which corresponds to a 
muon beam momentum of 90.622 GeV/c. This is the first spin resonance for muons. At this 
point, the fitting method loses sensitivity completely, since there will be no spin oscillations 
turn by turn. We now study the error £7/7 as a function of beam energy for P=-0.26 and 
PERR=0.5E-2 (keeping the magnetic field in the idealized storage ring to be 4.0 Tesla) 
as a function of muon beam energy that straddles the spin resonance. For initial muon 
collider physics, the interesting beam energies are 45.5 GeV (half the Z mass), 80.3 GeV 
(W threshold), 175 GeV (top threshold) as well as half the neutral Higgs mass, which could 
be as low as 55 GeV in some SUSY scenarios. We sample all electrons that have energies 
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TABLE II. Results of fits for ^7/7 as a function of polarization P and noise PERR. Also 
shown is the \ 2 of the fit for 1000 turns. 
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greater than half the muon energy. Figure |6] shows the variation of £7/7 as a function of 
muon beam energies that straddle these values. It can be seen that £7/7 first decreases as 
one gets close to the resonance and then blows up on the spin resonance. Figures (|?|-|TT|) 



show the fitted solutions superimposed on the simulated data for various momenta. Also 
shown side by side is the simulated data by itself. As one approaches the spin resonance, 
the oscillations slow down. It is nevertheless possible to fit the slowed down oscillations by 
a rapidly oscillating theoretical function to high accuracy on either side of the resonance. 
At the resonance, the oscillations die completely, which results in a large value of Sj/j. 
It may be possible to use this blow-up in £7/7 to find the spin resonance accurately and 
(paradoxically) determine 7 at resonance accurately. This would depend on the width of 
the spin resonance, an analysis of which would take us beyond the scope of this paper. 

B. Variation of ^7/7 as a function of beam energy spread 

We now calculate the variation of polarization as a function of turn number for an 
ensemble of muons with initial value of polarization P = -0.26 and values of momentum 
spread 5p/p varying from 0.02E-2 to 0.00125E-2. This variation is plotted in figure |13[ For 
the larger values of momentum spread, there is a significant degradation of polarization as a 
function of turn number, due to differential spin precession of the individual beam particles. 
We note that when the beam energy is at 175 GeV, the spin tune is significantly higher and 
the depolarization is more rapid. Despite this depolarization, there is enough information 
from the first few hundred turns to extract the excellent value of $7/7 for 175 GeV beam 
energy as shown in figure 0. 



Figure [14] shows the variation of the fractional energy resolution, £7/7 as a function 
of fractional beam energy spread for a muon beam with P = -0.26, with 41261 electrons 
sampled. There is little dependence of £7/7 on the momentum spread. This is due to the 
fact that the momentum spread is determined from the spin tune and not from the spin 
oscillation amplitude and the fact that the depolarization is not significant for the first few 
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hundred turns for any of the beam momentum spreads considered here. 



C. Optimization of the electron energy cut 

We now vary the cut on electron energy and study the dependence on £7/7 on the cut. 



Figure [15] shows the variation of £7/7 with the cut on individual electron energies for P 
= -0.26 for 41261 and 1146 electrons sampled. As shown in the Appendix, the fractional 
error on the average energy of electrons is much smaller than the fractional error on the 
total energy of electrons. It is possible to measure the average electron energy by counting 
the number of electrons going into the calorimeter with a scintillator array. However, the 
precession information is contained increasingly in the number of electrons rather than their 



average energy as we increase the electron energy cut. Figure |T| shows the variation of £7/7 
calculated from average as well as total electron energy as a function of the electron energy 
cut. For smaller values of the electron energy cut, the average method produces superior 
errors than the total energy method. However, with 40,000 electrons or more sampled a total 
energy method with a cut of 25 GeV or higher seems optimal. It should however be pointed 
out that the average energy method does not require a model for the rate of decay of muon 
intensity in the machine, which in practice could be a complicated function of turn number. 
As such the systematics associated with this would not be present in the average energy 



method. Figure |l6](a) shows the variation of the absolute value of C/F as a function of the 



electron energy cutoff for P = -0.26, where C and F are defined in equation |1.10| for both 
the total energy method and the average energy method. Figure [l6|(b) shows the fraction 
of electrons that lie above the electron energy cut as a function of the energy cut. The 
polarization for this sample is 0, since the electron energy fraction depends on polarization 



as well. Given the curves shown in figure |Lg, it should be possible to estimate the error in 
57/7 for a variety of conditions. 
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III. EFFECTS DUE TO DEPARTURES FROM THE IDEAL CASE 



So far we have considered a planar collider ring with uniform vertical magnetic field and 
no electric fields. The actual collider ring will depart from the ideal in three respects; a) It 
will have RF electric fields to keep the muons bunched, b) it will have radial horizontal 
magnetic fields experienced by partcles in an off-center trajectory at quadrupoles and at 
vertical correction dipoles, and c) it will have longitudinal magnetic fields due to solenoidal 
magnets in the interaction region(s). We now consider the effect due to each of these 
departures from the ideal. 

A. Electric fields 

Equation |1.2| implies that there is no spin precession due to longitudinal electric fields 
(/? x E — 0). RF electric fields are longitudinal, so there will be no precession due to the RF 
electric fields. At present there are no plans to install electrostatic separators to separate 
the beams. If and when this happens, one should consider the effect due to the transverse 
electric fields thus introduced. 

B. Effect of radial magnetic fields 

Particles which are off-axis at quadrupoles will experience radial as well as vertical mag- 
netic fields. Even though the net integral of these off-axis fields around the ring is zero, the 
spin rotation along a horizontal axis followed by spin rotation about a vertical axis (caused 
by a bend dipole) followed by a reverse rotation in the horizontal direction still produces a 
net effect since the rotations about the horizontal and vertical axes do not commute. The 
effects have been analyzed by Assmann and Koutchouk who show that this results in 
both a net spin tune shift < 5u > as well as a spread in tune o& v . 

< 8v >= C ^vl {n Q {Kl Q fal + n cv a* cv ) (3.1) 
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where v$ = 07 is the spin tune of the collider ring, uq are the number of quadrupoles 
with integrated gradient KIq, a y is the misalignment spread of the closed orbit at the 
quadrupoles, ncv is the number of vertical correction dipoles and aecv is the rms beand 
angle in the vertical correctors. The spread in tune is given by, 

< 5v > ,„ , 

<T8v = 3.2 

COSTTUq 



Table |T| shows the values for < 5i> > and ag u obtained by Assman and Koutchuk ||] for 



LEP. We compare this with to the current design for the 50 GeV muon collider ring ||. 
Including the low beta section, there are 70 quadrupoles with an RMS value of KIq = 0.27 
m _1 . The effects due to correction dipoles may be neglected in both the LEP and the muon 
collider cases. We assume a beam misalignment of 5mm at the quadrupoles, which is the 
same value used in the LEP calculation. This is probably being conservative. The tune shift 
for LEP corresponds to a shift in beam energy calibration of 3.0 KeV. The tune spread for 
LEP corresponds to a spread in beam energy calibration of 30 KeV. For the muon collider, 
the tune shift corresponds to a shift in beam energy calibration of -0.24 KeV and a spread 
of 1.46 KeV, both of which are negligible. The reason for the smallness of this effect for the 
muon collider is twofold. Since the circumference of the muon collider is smaller than LEP, 
there are fewer quadrupoles. Secondly, the muon is two hundred times more massive than 
the electron and has has a spin tune 07 that is smaller by the same factor. The spin tune 
shift depends on the the square of the spin tune. It should be noted that the above formulae 
are not valid for a fractional spin tune of 0.5. 

C. Solenoidal magnetic fields 

The experimental region will in all likelihood contain a solenoidal magnet. This solenoidal 
field, if uncorrected, will rotate the spin vector of the muons about the beam direction by a 



constant amount 9 S per turn, which can be derived using equation IO 



S = —(l + a)B s = -(l + a )^- (3.3) 
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where B s is the field due to the solenoid of length I, B is the dipole bending field of the ring 
of radius p. For a solenoid of 1.5 Tesla and length 6 meters, 9 S = 3.09 degrees for the planar 
storage ring parameters of table |I[ It can be shown analytically |J that this produces a spin 
tune shift bv given by 

v + 5u = —arccos ( cos (iris) cos (—) ) (3.4) 

7T \ 2 J 

yielding a spin tune shift 5u = -1.901E-5, or a fractional spin tune shift of bvjv = -3.45E-5. 
For a 50 GeV muon beam, this is a shift in energy calibration of -1.72 MeV. In LEP, a 
similar solenoid will have a much smaller fractional tune shift || , since the tune is 200 times 
larger for electrons. It is important to correct the effect due to the solenoids, since this is 
cumulative turn by turn. At LEP this is done by a series of vertical orbit correctors || 
followed by normal lattice followed by vertical orbit correctors of reverse polarity, which has 
the effect of rotating the spin by half the amount produced by the solenoid. A similar set of 
corrections is inserted after the solenoid to complete the correction. This method depends 
on a non-zero value of g — 2 and as such will be 200 times less effective for muons than 
for electrons, for any given magnet strength. The most effective method to correct for the 
solenoid is to surround it on either side by compensating solenoids of minimal radius large 
enough to allow the beam to go through. 



IV. CONCLUSIONS 

We have demonstrated that it is feasible to measure the energy of a 50 GeV muon 
collider to a few parts per million using the g — 2 spin precession technique, provided it is 
feasible to maintain a muon polarization of the order of P=0.25 in the ring for a thousand 
turns. In order to explore the Higgs resonance, it is necessary to measure the bunch by 
bunch variation in energy to a few parts per million. We have demonstrated that the g — 2 
technique is capable of doing so. It is still possible to tolerate a spin tune shift in the overall 
energy scale of a few percent, which will act only as a systematic error on the Higgs mass 
and width. 
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We would also like to note in passing that polarization information from a calorimeter 
of the type proposed here can be used in conjunction with a neutrino detector placed along 
the line of the neutrinos produced in association with the electrons to estimate the variation 
in the energy spectrum of the muon neutrinos and electron antineutrinos in the beam. 
Such information can be a valuable tool in untangling various possible neutrino oscillation 
scenarios. 

We intend to develop the method here by studying the propagation of polarized muons 
in a realistic 50 GeV collider lattice using the program COSY @ that takes into account 
non-linear effects in the dynamic aperture. Design and Monte Carlo studies will also be 
undertaken to develop the calorimeter detector needed. The authors would like to acknowl- 
edge useful conversations with Alain Blondel, Yaroslav Derbenev and Robert Rossmanith. 
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FIG. 4. (a)Energy detected in the calorimeter during the first 50 turns in a 50 GeV muon 
storage ring (points). An average value of P=-0.26 is assumed and a fractional fluctuation of 



0.5E-2 per point. The curve is the result of a MINUIT fit to the functional form in equation 1.1C. 

(b) The same fit, with the function being plotted only at integer turn values. A beat is evident. 

(c) Pulls as a function of turn number (d)Histogram of pulls. 
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FIG. 5. (a)Fractional error in £7/7 obtained from the oscillations as a function of polarization 
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FIG. 6. Fractional error in £7/7 obtained from the oscillations as a function of muon beam 
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FIG. 7. The figures on the left hand side show the simulated data with the fitted function 
superimposed for 50 turns. The figures on the right hand side show the simulated data and the 
fitted function at integer values of the turn number. The data shown are 60 GeV/c and 70 GeV/c 
muon momenta respectively. 
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FIG. 8. The figures on the left hand side show the simulated data with the fitted function 
superimposed for 50 turns. The figures on the right hand side show the simulated data and the 
fitted function at integer values of the turn number. The data shown are 80 GeV/c and 90 GeV/c 
muon momenta respectively. 
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FIG. 9. The figures on the left hand side show the simulated data with the fitted function 
superimposed for 50 turns. The figures on the right hand side show the simulated data and the 
fitted function at integer values of the turn number. The data shown are 90.622 GeV/c and 91.2 
GeV/c muon momenta respectively. The upper curve is on resonance. 
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FIG. 10. The figures on the left hand side show the simulated data with the fitted function 
superimposed for 50 turns. The figures on the right hand side show the simulated data and the 
fitted function at integer values of the turn number. The data shown are 100 GeV/c and 110 
GeV/c muon momenta respectively. 
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FIG. 11. The figures on the left hand side show the simulated data with the fitted function 
superimposed for 50 turns. The figures on the right hand side show the simulated data and the 
fitted function at integer values of the turn number. The data shown are 120 GeV/c and 130 
GeV/c muon momenta respectively. 



26 



x 10' 



Oscillations as a function of momentum 

x 10 3 



11/09/97 17.00 



5750 
5500 
5250 
5000 
4750 
4500 
4250 
4000 
3750 








175 GeV 




10 20 



30 40 50 
turn number 



10 20 



30 40 50 
Turn number 



FIG. 12. The figures on the left hand side show the simulated data with the fitted function 
superimposed for 50 turns. The figures on the right hand side show the simulated data and the fitted 
function at integer values of the turn number. The data shown are 175 GeV/c muon momentum. 
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FIG. 13. Variation of polarization as a function of turn number for 50 GeV muons with initial 
P =-0.26 and various values of Sp/p in an ideal collider ring. The bottom curve is for 175 GeV 
muons and shows a more rapid depolarization due to the higher spin tune. 
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FIG. 14. ^7/7 versus fractional beam energy spread for 50 GeV muons with PERR=.5E-2 
and P =-0.26 
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FIG. 15. The variation of £7/7 as a function of the electron energy cut for 41261 and 1146 
electrons P =-0.26. We fit the total energy in the calorimeter as well as the average energy per 
electron 
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FIG. 16. (a) The variation of C/F as a function of the electron energy cut for P =-0.26 for 
total energy method and average energy method, (b) The fraction of electrons that survive the 
energy cut as a function of the cut for P = 0. 
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V. APPENDIX 



A. Treatment of Errors 

We measure the total energy E of all electrons with individual energy e > 25 GeV in an 
electromagnetic calorimeter. Let N be the number of electrons sampled during a turn. N 
can fluctuate from sampling to sampling. Then 

TV 

E = Y i e i = N < e> (5.1) 
i=i 

<E> 2 <N> 2 <e> 2 <N> y <e> 2) V ' J 
where the variance cr 2 of the quantities e and E results from the kinematic distributions 
of those quantities and not from the measurement errors. The average of the individual 
electron energies is denoted by < e >. 

Let the calorimeter be such that it measures the true deposited energy E with a resolution 
e(E) such that 

where C,S and N represent the Constant, Sampling and Noise terms respectively. Let us 
assume that the measurement errors are Gaussian. Then, 

P(E m ) = J P(E)G(E, E m , e)dE (5.4) 

where E m is the measured energy and G(E,E m ,e) is a Gaussian of mean E and standard 
deviation e, which is a function of E and is written as e(E). From this it follows that the 
mean measured energy < E m > and the mean squared measured energy < E^ > are given 
by 

<E m >=J E m P(E m )dE m = J E m dE m J P(E)G{E, E m , e)dE 

(5.5) 

= / P(E)dE J E m G(E, E m , e)dE m = J P(E)dE x E =< E > 



The equation |5.5| states that the mean value of any distribution given by P(E) is the same 



as that of the smeared distribution P(E m ) provided the smearing function is such that 
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the average of the smeared values for any given true value E equals the true value ( a 
property satisfied by Gaussians) and the integration is carried over the full range of the 
variables. As an aside, in High Energy Physics, we measure steeply falling spectra that are 
smeared by measurement errors. Provided there is no arbitrary lower cut-off in the measured 
spectra (such as a trigger threshold), the above result would be valid, even for non-Gaussian 
resolutions. For the muon collider, the cut off in selected electrons of 25 GeV is imposed by 
momentum selection that is independent of the calorimetry. So the above result would still 
be valid. Similarly, one can compute < E^ > 

< E 2 m >= JE 2 m P(E m )dE m = JE 2 m dE m JP(E)G(E,E m ,e)dE 

(5.6) 

= J P(E)dE J ElG(E,E m ,e)dE m = J P(E)dE x (E 2 + e 2 ) =< E 2 > + J P(E)e 2 (E)dE 
From this it follows that the variance of the measured energy a% is given by 

<>k. =(T e + J P(E)e 2 (E)dE » o% + e 2 (< E >) (5.7) 

where the last approximation results from assigning the average measurement resolution to 
the resolution at the average energy. This then leads to 

°i °l e 2 (<E>) 



<E m > 2 <E> 2 < E >< 



Using equation f[2\ and 5.3 leads to 



a\ j_ a 2 e ^ 2 S 2 Af 2 

<E^> 2 ^N {1 + T7^ )+ + N < e > + N 2 < e > 2 (5 ' 9) 

From the above equation, it is obvious that the calorimeter must be such that the constant 

term C must be negligible for the fractional resolution to scale inversely with the number iV 

of electrons collected. The noise term can be neglected for large enough iV since it goes as 

iV -2 . With these assumptions, one gets 

al - (5.io) 



<E m > 2 N < e > 2 <e> 
For a 50 GeV muon beam, the values of < e > and a e are 34.05 GeV and 6.046 GeV 
respectively for electrons with e > 25 GeV. The ratio, a e j < e > is to a good approximation 
independent of muon energy. This then leads to the following error formula. 
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ol 1 , S 2 

— (1 + 0.03153 + ——) (5.11) 



<E m > 2 N x ' 34.05 
Sampling terms of 0.15 GeV 1 / 2 or better are easy to obtain in electromagnetic calorimeters. 



This leads to 

al 1 



< E m > 2 N 



1 + 0.03153 + 0.000661) (5.12) 



i.e. the sampling term can be neglected when compared to the fluctuation in the true 
electron energies. So if the fractional measurement error is PERR = ^. m > is specified, the 
equivalent number of electrons is given by 

1-03153 ,ri^ 

NK Tmm (5 - 13) 

In other words, PERR=.5E-2,1.0E-2,2.0E-2 and 3.0E-2 implies 41261, 10315, 2579, and 
1146 electrons sampled. If in practice we sample 100,000 electrons, this leads to a value of 
PERR=0.3212E-2. In order for this good a resolution to be meaningful, the constant term 
C has to be below this order of magnitude. 



B. Using averages 

Equation 1572] holds for the total energy E in the calorimeter. If however, one also 



measures the total number of particles entering the calorimeter (using a scintillator system 
for example, that counts minimum ionizing particles), then for each turn one can measure 
the average energy E av of electrons. The fractional error on E av does not contain a term 
due to the fluctuation of the number of electrons entering the calorimeter, being given by 



(5.14) 



< E av > 2 <e> 2 <N> < e > 2 

with < E av >=< e >. For a fractional error of PERR in E av , the equivalent number of 
electrons sampled would be given by 

N **TPEim (5 - 15) 
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With this method, PERR=.5E-2,1.0E-2,2.0E-2 and 3.0E-2 implies 1261, 315, 79, and 35 
electrons sampled, assuming no error in the measurement of N. If we sample 100,000 elec- 
trons, the fractional error in the average would be 0.561E-3. For this error to be meaningful, 
the sampling term would have to be of this order of magnitude. 
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Machine 


Spin tune vq 


Quadrupoles 


RMS Kl Q 
meters -1 


Gy 

meters 


8v 




46 GeV LEP 


100.47 


» 600 


0.032 


0.5E-3 


5.7E-6 
= 3KeV 


6.1E-5 
= 30KeV 


50 GeV Muon Collider 


0.5517 


70 


0.274 


0.5E-3 


-0.26E-8 
= -0.24KeV 


1.66E-8 
= 1.46KeV 



TABLE III. Predictions for spin tune shift bv and spread in spin tune shift o& v caused by 
quadrupoles for LEP compared to the 50 GeV muon collider ring 
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